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ABSTRACT 

We present three different optical and near-infrared (NIR) data sets for 
evolved stars in the Galactic Globular Cluster u> Centauri . The comparison 
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between observations and homogeneous sets of stellar isochrones and Zero- Age 
Horizontal Branches provides two reasonable fits. Both of them suggest that 
the so-called anomalous branch has a metal-intermediate chemical composition 
(—1.1 < [Fe/H] < —0.8) and is located ~ 500 pc beyond the bulk of uj Cen stars. 
These findings are mainly supported by the shape of the subgiant branch in four 
different color-magnitude diagrams (CMDs). The most plausible fit requires a 
higher reddening, E(B — V) = 0.155 vs. 0.12, and suggests that the anomalous 
branch is coeval, within empirical and theoretical uncertainties, to the bulk of 
uj Cen stellar populations. This result is supported by the identification of a 
sample of faint horizontal branch stars that might be connected with the anoma- 
lous branch. Circumstantial empirical evidence seems to suggest that the stars 
in this branch form a clump of stars located beyond the cluster. 

Subject headings: globular clusters: general - - globular clusters: individual 
(uj Centauri) 

1. Introduction 

The peculiar Galactic Globular Cluster (GGC) uj Centauri (NGC5139) is currently 
subject to substantial observational efforts covering the whole wavelength spectrum. This 
gigantic star cluster, the most massive known in our Galaxy, has (at least) three separate 
stellar populations with a large undisputed spread in age, metallicity (Fe, Ca) and kinematics 
(e.g. Hilker & Richtler 2000; Ferraro et al. 2002; Smith 2004). According to recent 
abundance measurements based on 400 medium resolution spectra collected by Hilker et al. 
(2004), it seems that the spread in metallicity among the uj Cen stars is of the order of 
one dex (—2 < [Fe/H] < —1). The metallicity distribution shows three well-defined peaks 
around [Fe/H] = —1.7, —1.5 and —1.2 together with a few metal-rich stars at ~ —0.8. On 
the basis of high-resolution spectra, it has been suggested by Pancino (2004) that stars in 
the anomalous branch (Lee et al. 1999) might be more metal-rich, with a mean metallicity 
[Fe/H] ~ —0.5. A similar metal-rich tail was also detected by Norris et al. (1996) and 
by Suntzeff & Kraft (1996, and references therein). In the absence of conformity in the 
literature for the uj Cen RGB names, we here refer to the metal-poor component as uj 1, to 
the metal-intermediate as uj2, and to the anomalous branch as uj3. 

Results from new observations seem to pose as many new questions as they answer. It 
has been suggested (e.g. Lee et al. 1999; Hughes et al. 2004) that the uj3 population might 
be significantly younger than uj2. However, in a recent investigation based on Very Large 
Telescope (VLT) and Hubble Space Telescope (HST) data, Ferraro et al. (2004) found that 
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the lo 3 population is at least as old as the uj 2 population and probably a few Gyrs older. 
The observational scenario was further complicated by the results brought forward by Bedin 
et al. (2004) on the basis of multi-band HST data. They not only confirmed a bifurcation 
along the Main Sequence, but also found a series of different Turn-Offs (TOs) and sub- 
giant branches (SGBs). The proposed explanations in the literature for these peculiar stellar 
populations in uj Cen are many and include, e.g., increased He content, or a separate group 
of stars located at a larger distance (Ferraro et al. 2004; Bedin et al. 2004; Norris 2004). 
The latter hypothesis could be further supported by the occurrence of a tidal tail in uj Cen 
but we still lack a firm empirical detection (see Leon et al. 2000; Law et al. 2003). The 
formation history and composition of uj Cen thus form a complex puzzle that is being slowly 
pieced together by investigations based on the latest generation of telescopes. A well-known 
promising technique is to study evolved stars simultaneously in optical and NIR bands to 
limit subtle errors due to the absolute calibration, to crowding in the innermost regions, and 
to reddening corrections. The main aim of this Letter is to investigate whether different 
assumptions concerning the spread in age and in chemical composition, or differences in 
distance and in reddening, account for observed optical (BRI) and NIR ( JK) CMDs. 



2. Observations and data reduction 



Near-infrared J, K s images of uj Cen were collected in 2003 with SOFI at the New 
Technology Telescope of ESO, La Silla. The seeing conditions were good and range from / /6 
to l'/l. Together with additional data from 2001, available in the ESO archive, we end up with 
a total NIR sample of 92 J- and 135 i^ s -band images that cover a 14 x 14 arcmin 2 area centred 
on the cluster. These data were reduced with DA0PH0TII/ALLFRAME following the same 
technique adopted by Dall'Ora et al. (2004). The NIR catalogue includes ~1.4 x 10 5 stars. 
With FORS1 (standard-resolution mode) on the ANTU/VLT telescope, we also collected 
optical UVI images in 1999. These data are from a 2 x 2 pointing mosaic centred on the 
cluster that covers an area slightly larger than 13 x 13 arcmin 2 . Exposure times range from 
10 (J) to 30 s (U) and seeing conditions were better than l'.'O. These data have also been 
reduced using DA0PH0TII/ALLFRAME and the photometric catalogue includes ~5 x 10 5 stars. 
Optical F435W and F625W (hereinafter B and R bands) data were retrieved from the HST 
archive. These images were collected with the ACS instrument in 9 telescope pointings, each 
of which provides a 5i?-image pair with exposure times of 12 (B) and 8 s (R). The field 
covered by these data is ~9 x 9 arcmin 2 , centred on the cluster. The ACS data were reduced 
with ROMAFOTwo and provide BR photometry for ~ 4 x 10 5 stars. The photometry was 
kept in the Vega system (see e.g., http://www.stsci.edu/hst/acs/documents). A detailed 
description of observations and data reduction will be given in a future paper. Here we only 
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wish to mention that raw frames were prereduced using standard IRAF procedures and, in 
addition, the FORS1 images were corrected for amplifier cross-talk by following Freyhammer 
et al. (2001). To improve the photometric accuracy, carefully chosen selection criteria were 
applied to pinpoint a large number (>100) of point-spread function stars across the individual 
frames, and several different reduction strategies were used to perform the photometry over 
the entire data set. The observed fields were combined to the same geometrical system using 
iraf . immatch and DAOMATCH/DAOMASTER. The absolute photometric calibration of ground- 
based instrumental magnitudes was performed using standard stars observed during the 
same nights. The typical accuracy is 0.02 -0.03 mag for both optical and NIR data. 



3. Results and Discussion 

Figure 1 shows CMDs based on the five different photometric bands. The plotted stars 
were selected from individual catalogues by using the 'separation index' sep introduced by 
Stetson et al. (2003), since crowding errors dominate the photometric errors. The adopted 
sep ranges from 6.5 (NIR data) to 8 (ACS data), which corresponds to stars having less 
than 0.3% of their measured light contaminated by neighbour stars; the higher the sep for 
a star, the less the severity of the crowding by neighbours. Error bars plotted at the left 
side of each panel account for photometric and calibration errors in magnitude and color. 
The input physics adopted in our evolutionary code has been discussed in detail in a series 
of papers (Cariulo et al. 2004; Cassisi et al. 1998). Here, we point out that the stellar 
models (partly available at http://astro.df.unipi.it/SAA/PEL/ZO.html) account for atomic 
diffusion, including the effects of gravitational settling, and thermal diffusion with diffusion 
coefficients given by Thoul et al. (1994). The amount of original He is based on a primordial 
He abundance Yp = 0.23 and a He-to-metals enrichment ratio of AY/AZ ~ 2.5 (Pagel & 
Portinari 1998; Castellani & DeglTnnocenti 1999). We adopt the solar mixture provided by 
Noels & Grevesse (1993). For details on the calibration of the mixing length parameter and 
on model validation, we refer to Cariulo et al. (2004) and Castellani et al. (2003). To avoid 
deceptive uncertainties in the comparison between theory and observations, the predictions 
were transformed into both the BR Vega system and the UK Johnson- Cousins bands by 
adopting the atmosphere models provided by Castelli et al. (1997) and Castelli (1999). 

Data plotted in Fig. 1 display the comparison between observations and a set of four 
isochrones (solid lines) with the same age (12 Gyr) and different chemical compositions (see 
labels). A mean reddening of E(B - V) = 0.12 (E(B - V) = 0.11 ± 0.02, Lub (2002)) was 
adopted and a true distance modulus of /i = 13.7 ((m — M)y = 14.05 ± 0.11; Thompson 
et al. (2001), (m — M)k = 13.68 ± 0.07, Del Principe et al. 2004, private communication). 
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Extinction parameters for both optical and NIR bands have been estimated using the extinc- 
tion model of Cardelli et al. (1989). The metal-poor and the metal- intermediate isochrones 
supply, within current empirical and theoretical uncertainties, a good fit to the bulk of RGB 
and HB stars. Close to the RGB tip, the isochrones are slightly brighter than the observed 
stars, which is caused by the adopted mixing-length parameter (a = 2). The SGB and the 
lower RGB are only marginally affected by this parameter, since the empirical isochrone 
calibration is based on these evolutionary phases (Cariulo et al. 2004). However, the most 
metal-rich isochrone appears to be systematically redder than the stars at the base of the 
uj 3 branch 1 (green dots) and fainter than the SGB stars. Clearly, the shape of the uj 3 SGB 
does not support metal-abundances > Z = 0.004, since TO stars for more metal-rich popu- 
lations become brighter than SGB stars. To further constrain this evidence, Fig. 2 shows the 
comparison between observed HB stars (the entire sample includes more than 2,300 objects) 
and predicted Zero Age Horizontal Branches (ZAHBs) with a progenitor age of 12 Gyr and 
different chemical compositions. For the sake of the comparison, the objects located between 
hot HB stars and RGB stars were selected in the B — R, B plane and plotted as blue dots. 
Theoretical predictions plotted in this figure show that the more metal-rich ZAHB appears 
to be systematically brighter than the observed HB stars for B — J B — K & 0.5. More- 
over, the same ZAHB crosses the RGB region, thus suggesting that the occurrence of HB 
stars more metal-rich than Z=0.002 should appear as an anomalous bump along the RGB. 
Data plotted in these figures disclose that the stellar populations in uj Cen cannot be ex- 
plained with a ranking in metal abundance. In line with Ferraro et al. (2004) and Bedin et al. 
(2004), we find that plausible changes in the cluster distance, reddening, age, and chemical 
composition do not supply a reasonable simultaneous fit of the uj 1, uj2, and uj3 branches. 

To further investigate the nature of the a; 3 branch, we performed a series of tests by 
changing the metallicity, the cluster age, and the distance. A good fit to the anomalous 
branch is possible by adopting the same reddening as in Fig. 1, a true distance modulus of 
fi = 13.9, and an isochrone of 15.5 Gyr with Z = 0.0025 and Y = 0.248. Figure 3 shows that 
these assumptions supply a good fit in both the optical and the NIR bands, and indeed the 
current isochrone properly fits the width in color of the sub-giant branch and the shape of a 
good fraction of the RGB. Moreover, the ZAHB for the same chemical composition (dashed 
line) agrees quite well with the faint component of HB stars. Note that more metal-poor 
ZAHBs at the canonical distance do not supply a reasonable fit of the yellow spur stars with 
15 < B < 15.5. The yellow spur is visible in all planes, and in the NIR (c,d) it even splits 
up in brightness due to a stronger sensitivity to effective temperature. The identification 



1 Stars of the w3 branch have been selected following the original detection in the B — R,R plane by 
Ferraro et al. (2004, see their Fig. 2). 
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of the entire sample was checked on individual images and, once confirmed by independent 
measurements, the data indicate a separate HB sequence for the u 3 population. 

Although the u3 fit may appear good, it implies an increase in distance of A/i = 0^2 
and an unreasonable ~ 4 Gyr increase in age. This estimate is at variance with the absolute 
age estimates of GGCs (Gratton et al. 2003) and with CMB measurements by W-MAP 
(Bennett et al. 2003). The discrepancy becomes even larger if we account for the fact that 
this isochrone was constructed by adopting a He abundance slightly higher (Y = 0.248, vs 
0.238) than estimated from a He-to-metals enrichment ratio AY/AZ = 2.5. This increase 
implies a decrease in age of ~ 1 Gyr. Moreover, we are performing a differential age estimate, 
and therefore if we account for uncertainties: in the input physics of the evolutionary models 
(e.g. equation of state, opacity); in the efficiency of macroscopic mechanisms (like diffusion); 
in model atmospheres applied in the transformation of the models into the observational 
plane; and in the extinction models, we end up with an uncertainty of ~ 1 — 2 Gyr in cluster 
age (Castellani & Degl'Innocenti 1999; Krauss & Chaboyer 2003). 

Owing to the wide range of chemical compositions and stellar ages adopted in the 
literature for explaining the morphology of the to 3 branch, we decided to investigate whether 
different combinations of assumed values of distance, chemical composition, and reddening 
may also simultaneously account for the u 3 branch and the HB stars. We found that two 
isochrones of 13 Gyr for Z = 0.0015 and Z = 0.003 bracket the a; 3 stars (see Fig. 4), within 
empirical and theoretical uncertainties. The fit was obtained using the same true distance 
modulus adopted in Fig. 3, together with a mild increase in reddening, E(B — V) = 0.155. 
Once again theory agrees reasonably well with observations in all color planes. Moreover, 
data plotted in Fig. 5 show that the predicted ZAHB with Z=0.0015 and Z=0.003, for the 
adopted distance modulus and reddening correction, account for the yellow spur stars. The 
same figure shows a sample of 53 RR Lyrae stars selected from the variable-star catalogue 
by Kaluzny et al. (2004) for which we have a good coverage of J- and .fT-band light curves 
(Del Principe et al. 2004, private communication). The RR Lyrae stars only account for a 
tiny fraction of the yellow-spur stars. This finding together with the detection a well-defined 
sequence in panels b,c, and d indicate that this spur might be the HB associated with the 
ui 3 population. Note that a few of these ' lj 3-HB' spur stars have also been detected by Rey 
et al. (2004, see their Fig. 7, for V « 14.75) and by Sollima et al. (2004a, see their Figs. 6 
and 12). 
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4. Final remarks 

We have presented a new set of multi-band photometric data for the GGC uj Cen 
and — in agreement with previous findings in the literature — we find no acceptable fit to 
the different stellar populations for a single distance, reddening, and age. We found two 
reasonable fits for the a>3 stars: (1) by adopting a 0.2 higher distance modulus (~ 500 pc), 
a metal-intermediate composition [Z = 0.0025, Fe/H « —0.9), and an unreasonable increase 
in age of ~4 Gyr; or (2) for the same A/i = 0.2 shift, an increase in the reddening, metal- 
intermediate chemical compositions (0.0015 < Z < 0.003, —1.1 < Fe/H < —0.8), and an 
age that, within current uncertainties, is coeval with the bulk of the uj Cen stars. Current 
findings indicate that uj 3 stars are not significantly more metal-rich than Z = 0.003. This 
evidence is supported by the shape of the to 3 SGB, as already suggested by Ferraro et al. 
(2004) and by the fit of HB stars. We are in favour of the latter solution for the following 
reasons: 

• The difference in distance between the a; 3 branch and the bulk of uj Cen stars is 
of the order of 10%. This estimate is 3-4 times smaller than the estimate by Bedin et al. 
(2004) and in very good agreement with the distance of the density maxima detected by 
Odenkirchen et al. (2003) along the tidal tails of the GGC Pal 5. Moreover, recent N-body 
simulations (Capuzzo Dolcetta et al. 2004) indicate that clumps along the tidal tails can 
approximately include 10% of the cluster mass. 

• We found that by artificially shifting the uj 3-branch stars to account for the assumed 
difference in distance and reddening, they overlap with the uj2 population. It has been 
recently suggested by Piotto et al. (2004), on the basis of low-resolution spectra, that the 
bluer main sequence detected by Bedin et al. (2004) is more metal-rich than the red main 
sequence. Unfortunately, current ACS photometry is only based on shallow ACS exposures, 
and therefore we cannot properly identify in our data these stellar populations located in the 
lower main sequence. The same outcome applies to the suspected extremely-hot HB progeny 
of the bluer main sequence, since they have not been detected in the NIR bands. 

• Current preliminary findings support recent N-body simulations by Chiba & Mizutani 
(2004) and by Ideta & Makino (2004). In particular, the latter authors found, by assuming 
that the progenitor of Omega Centauri is a dwarf galaxy, that more than 90% of its stellar 
content was lost during the first few pericenter passages (see their Fig. 2). 

• uj Cen reddening estimates in the literature cluster around E(B — V) = 0.12 ± 0.02. 
However, the map from (Schlegel et al. 1998) indicates reddening variations of 0.02 across 
the body of the cluster while, more importantly, 2MASS data (Law et al. 2003) show a 
very clumpy reddening distribution at distances beyond 1° (100 pc) from the cluster centre, 
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with large variations AE(B — V) — 0.18 across a 4 degrees 2 area. However, the inference 
of reddening between the main body of uj Cen and the supposed background population is 
very surprising, uj Cen lies at the comparatively low Galactic latitude of +15°, and the 
reddening variations seen in the Schlegel and 2MASS maps likely originate in the foreground 
interstellar material of the Galactic disk. Any interstellar material behind uj Cen must lie 
at least 1.4 kpc from the Galactic plane, and therefore would most likely be associated with 
uj Cen itself. Smith et al. (1990) have reported a significant detection of H I in the direction 
of uj Cen , blueshifted by ~ 40 km s -1 with respect to the cluster velocity. However, their 
interpretation is that this gas is associated with the northern extension of the Magellanic 
Stream far beyond the cluster. It would be a remarkable coincidence if this interstellar 
material happened to lie between the cluster and a clump in the tidal tail at a 10% greater 
distance than the main cluster body, while traveling at the quoted relative speed. However, 
it is worth noting that Smith et al. (1990) estimated for this cloud a column density of 
Nh ~ 3 x 10 18 atmos/cm 2 , and therefore a reddening E(B — V) ~ 0.07 (Predehl & Schmitt 
1995) that is at least a factor of two larger than required by our best fits. This point is 
crucial for the proposed explanation and needs to be further investigated. 

Finally, we note that the comparison between predicted ZAHBs and HB stars indicates 
that the occurrence of an old stellar population with Z > 0.002 would imply the occurrence 
of an anomalous clump along the RGB. In fact, more metal- rich, red HB stars cover the 
same color range of metal-poor RGB stars. The detection of such a feature along the RGB 
can supply robust constraints on the progeny of the uj3 stellar population. These results, 
when independently confirmed, would suggest that the uj 3 branch might be a clump of stars 
located 500 pc beyond the bulk of the cluster. No firm conclusion can be drawn on the basis 
of current data, although this evidence together with the increase in radial velocity among 
uj3 stars measured by Sollima et al. (2004b) and numerical simulations recently provided 
by Capuzzo Dolcetta et al. (2004) indicates that it could be a tidal tail. 
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Fig. 1. — Optical (panels a,b), and NIR (panels c,d) CMDs for selected subsamples of the 
detected stars, compared with a set of 12 Gyr isochrones (solid lines) at different chemical 
compositions (see color coding). The adopted true distance modulus and cluster reddening 
are \x = (m — M) = 13.7 and E(B — V) = 0.12, respectively. The number of stars selected 
(NS) is also indicated. 
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Fig. 2. — Same as in Fig. 1, but the comparison between theory and observations is focused 
on Horizontal Branch stars. Blue objects mark stars located between hot HB stars and RGB 
stars. They have been selected in B — R,B plane. 



-14- 




-10 T 2 3 4 5 6 
B-R/B-l/B-J/B-K [mag] 



Fig. 3. — Same as in Fig. 1, but compared to a 15.5 Gyr isochrone for the on 3 population, 
constructed by adopting Z = 0.0025, and Y = 0.248. The fit was performed by adopting 
a true distance \i = 13.9, and a reddening correction EiB — V) = 0.12. Note that the 
corresponding ZAHB matches the selected HB stars (blue objects). 
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Fig. 4. — Same as Fig. 3, but compared to 13-Gyr isochrones constructed by adopting dif- 
ferent chemical compositions (see labels) and a higher reddening. 
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Fig. 5. — Same as in Fig. 4, but the comparison between theory and observations is focused 
on Horizontal Branch stars. Blue objects mark stars located between hot HB stars and RGB 
stars, while red dots (panels c,d) in the HB region display RRLyrae stars for which we have 
accurate mean NIR magnitudes. 



